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Abstract 

O ■ We prove an analog of Levinson's theorem for scattering on a weighted (to + 1)- vertex graph 

with a semi-infinite path attached to one of its vertices. In particular, we show that the number 
of bound states in such a scattering problem is equal to to minus half the winding number of 
the phase of the reflection coefficient (where each so-called half-bound state is counted as half a 
bound state). 
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^ ■ 1 Introduction 
O 

t^- | Levinson's theorem [7] , a classic result in quantum scattering theory, reveals a surprising connection 

between the scattering and bound states of a potential on a half line: the bound states of the 
potential are counted by the winding of the phase of the reflection coefficient. 

Here we consider the setting of scattering on graphs. In this model, scattering occurs on a 
discrete object, namely a (weighted) finite graph. Semi-infinite paths are attached to some of the 
vertices, allowing for incoming and outgoing wave packets. This setting was introduced (in the case 
of two semi- infinite paths) as a framework for developing algorithms for quantum computers [4]. 
Subsequently, scattering on graphs was used to discover a quantum algorithm for evaluating game 
trees [3], resolving a longstanding open question in quantum query complexity. The generalization 
to many semi-infinite paths was used to design a model of universal quantum computation [2]. 

In fact, discrete analogs of quantum scattering theory were considered prior to the advent of 
quantum computation. In the most natural discrete version of a potential on a half line, scattering 
occurs on a semi-infinite path of vertices with weighted self-loops representing the potential. Case 
and Kac proved an analog of Levinson's theorem for this setting [1], and Hinton, Klaus, and Shaw 
generalized their result to apply to systems with bound states of a special kind, called half-bound 
states [6]. 

In this paper, we prove an analog of Levinson's theorem for scattering on an arbitrary weighted 
finite graph with a semi-infinite path attached to one of its vertices. In particular, if the original 
finite graph contains m vertices, we show that the number of bound states is equal to m minus 
half the winding number of the phase of the reflection coefficient, where each half-bound state is 
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counted as half a bound state. The main idea of the proof is to consider the analytic continuation 
of the reflection coefficient from the unit circle to the rest of the complex plane. Using the Cauchy 
argument principle, we relate the winding of the phase of the reflection coefficient around the unit 
circle to the zeros and poles of the analytically continued reflection coefficient inside the unit circle. 
Finally, we relate these zeros and poles to bound states. 

The remainder of the article is organized as follows. In Section 2 we introduce the model of 
scattering on graphs and describe scattering states and bound states. In Section 3 we prove the 
main result. We conclude in Section 4 with a discussion of possible directions for future work. 



2 Scattering on graphs 

Consider an infinite path of vertices, each corresponding to a basis state \x) for x € Z, where vertex 
x is connected to vertices xdbl. The eigenstates of the adjacency matrix of this graph, parametrized 
by k 6 [— vr,7r), are the momentum states \k) with 

(x\k) = e ikx (1) 

for each x £ Z. These states are normalized so that (k\k') = 2ir 5(k — k'). Letting -ff pa th denote the 
adjacency matrix of this infinite path, we have 

{x\H path \k) = (x - l\k) + (x + l\k) (2) 
= 2cos k(x\k), (3) 

which shows that \k) is an eigenstate of -ff pat h with eigenvalue 2 cos A:. 

Now let G be an (m + l)-vertex graph, and create an infinite graph by attaching a semi-infinite 
path to the first vertex of G. We refer to the remaining m vertices of G as internal vertices. Label 
the basis states for vertices on the semi-infinite path as \x), where x = 1 for the vertex in the 
original graph and x = 2, 3, . . . moving out along the path. 

In general, we may assign complex weights to the directed edges of G, subject to the Hermitic- 
ity constraint that edges in opposite directions have conjugate weights. We denote the weighted 
adjacency matrix of G by the block matrix 



a &t' 
b D 



(4) 



where a £ R, b € C m (a column vector), and D = € C' mxm . A semi-infinite path is attached to 
the first vertex of G, with all edges on the semi- infinite path having weight 1. The adjacency matrix 
of the resulting infinite weighted graph is denoted H. See Figure 1 for a graphical representation 
of this Hamiltonian. Note that the setting considered in Refs. [1,6], restricted to the case where 
the potential is only nonzero over a finite interval, is the special case where &t = (1,0, ... ,0) and 
D is a tridiagonal matrix with off-diagonal entries equal to 1. 

To analyze scattering on G, it is useful to determine the eigenstates of H. Next we describe the 
two families of eigenstates, scattering states and bound states. 

2.1 Scattering states 

For each \ 7rZ, there is a scattering state |sc(fe)) of momentum k. This state has the form 

(x|sc(Jfe)> = e~ [kx + R(e [k ) e ikx (5) 
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Figure 1: The scattering Hamiltonian. 



on the semi-infinite path, where R(e lk ) is called the reflection coefficient. We can view R as a 
function of e lk since it is 27r-periodic in k. By a similar calculation as in (3), such states satisfy the 
eigenvalue condition 

H\sc(k)) = 2cosk\sc(k)) (6) 

at vertices of the semi-infinite path. Requiring (6) to hold at the m + 1 vertices of G gives a system 
of m + 1 linear equations that determines the reflection coefficient R(e lk ) and the amplitudes of 
|sc(fc)) at the m internal vertices of G. 

Given the Hamiltonian H defined in terms of the weighted adjacency matrix (4), it is straight- 
forward to solve for the reflection coefficient. The eigenvalue conditions for the scattering states 
|sc(fc)) at the vertices of G can be collected into the vector equation 

a &t\ / e " ifc + R(e ik ) e ik \ fe~ 2ik + R{e ik ) e 2ik \ _ , (e~ ik + R(e ik ) e ik \ 
b Dj{ V(e ifc ) ) + { ) ~ 2cosk { ^( e *) ) ' W 

where ij){e lk ) € C m is the vector of amplitudes at the m internal vertices of G in the state |sc(/c)). 
The lower block gives 

^(e ik ) = (e" ifc + R(e ik ) e ik )(2 cos k - D)~ l b. (8) 
Applying this identity to the upper block of (7) and solving for R(e lk ), we find 1 

where 



Q(e ik ) :=l-e ik (a + C(e ik )) (10) 
C(e ik ) :=6 t (2cosA;- J D)- 1 6. (11) 

It is easy to see that |ii(e lfc )| = 1 since Q(e lk )* = Q(e~ lk ). 

x Note that the expression (9) depends on the choice (5) for the form of the scattering states. In particular, it 
assumes a particular convention for the overall phases of these states. While the stated form will turn out to be 
convenient for our purposes, other natural choices can be obtained by taking states of the same form, but numbering 
vertices on the semi- infinite path starting from an integer other than 1. Such a choice modifies the reflection coefficient 
by factors of e lfc , which ultimately modifies the statement of Levinson's theorem. 



3 



The definition of R{e lk ) as a parameter describing scattering states only makes sense for fcgl. 

Nevertheless, it is convenient to extend the definition of R to the whole complex plane. For all 
z G C we let 

^ ); =-W (12) 

where 

Q(z) :=1- z(a + C(z)) (13) 

C{z) :=b\z + z- 1 - D)~ l b. (14) 

With z = e lk , these definitions agree with (9)-(ll). In fact, R(z) is the analytic continuation of 
R(e lk ) from the unit circle r := {z € C: \z\ = 1} to C. Note that for z ^ T, we may have |-R(z)| ^ 1. 

2.2 Bound states 

Scattering states alone may not span the entire state space. In general, to obtain a complete basis, 
we must also include bound states. 

For k > 0, a state |bd ± (K)) of the form 

(x|bd ± ( K )) = (±e- K f (15) 

on the semi-infinite path satisfies 

i^|bd ± (K)) = ±2 cosh k |bd ± ( K )) (16) 

at vertices of the semi-infinite path. If some choice of the amplitudes at the m internal vertices 
of the graph also satisfies the condition (16) on the vertices of G, then |bd ± («)) is an eigenstate 
of H with eigenvalue ±2 cosh k. The m eigenvalue conditions at the internal vertices of G give 
linear equations for the m amplitudes at internal vertices in terms of k. The remaining eigenvalue 
condition at the first vertex of G gives a transcendental equation in k with a discrete set of solutions, 
each corresponding to a bound state. 

The states |bd ± (K)) can be characterized as follows. Such a state satisfies the vector equation 



&t\ f ±e~ K \ fe- 2K \ , „ , /±e 



b Djy^( K )j + { o ) =±2coshK {^( K )J (17) 

where ^(k) € C m is a vector of amplitudes at the m internal vertices of G. The lower block gives 

^±( K ) = e~ K (2cosh«;=F J D)" 1 6. (18) 

Using this in the upper block of (17) gives Q(±e~ K ) = 0, where Q(z) is given in (13). In other 
words, bound states correspond to values of k for which ±e _K is a root of Q{z). 

The states |bd =t («;)) with k > have amplitudes that decay exponentially along the semi- 
infinite path, so we refer to them collectively as evanescent bound states. However, in the preceding 
discussion, we have allowed the possibility that k = 0. States of the form |bd =t (0)) may or may not 
exist depending on whether Q(±l) = 1 =F (a + C(±l)) = 0. Using the terminology of Ref. [6], we 
refer to such states as half-bound states. This reflects that they are, in a sense, halfway between 
scattering and bound states: when they exist, we have |bd + (0)) = |sc(0)) and |bd~(0)) = |sc(7r)). 
It also reflects the role they play in Levinson's theorem. 
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Finally, there may be bound states with zero amplitude on the semi-infinite path. We refer to 
such states as confined bound states (they have also been called bound states of the second kind in 
the terminology of Ref. [8]). A state |cd) with zero amplitude on the semi- infinite path satisfies 
the eigenvalue equation H\cd) = A|cd) on the vertices of the semi- infinite path for any value of A. 
To satisfy the eigenvalue equation at the internal vertices of G, a state |cd) must satisfy the vector 
equation 



The lower block says that ip is an eigenvector of D with eigenvalue A, and the upper block says 
that b^ip = 0. Thus, confined bound states correspond to eigenvectors ip of D that also satisfy the 
further condition of orthogonality to b. Note that if A is a degenerate eigenvalue of D, then there 
may be multiple confined bound states corresponding to A. In general, the number of confined 
bound states corresponding to A is the dimension of the subspace {ip: Dip = Xip, b^ip = 0}. If b is 
orthogonal to the A-eigenspace of D, this dimension is simply the multiplicity of A; otherwise it is 
the multiplicity of A minus 1. 

It can be shown that the scattering states {|sc(A;)} : k € (— vr, 0)}, together with the bound states 
{|bd ± (K)) : k > 0, Q(±e~ K ) = 0} and confined bound states {|cd)}, form a complete basis for the 
state space [5]. 

3 Levinson's theorem 

In order to state Levinson's theorem, we first introduce the concept of the winding number of a 
complex-valued function around a curve in the complex plane. Let 7 be a simple, closed, positively- 
oriented (i.e., counterclockwise) curve in C, and let /: C — > C be a meromorphic function with 
no zeros or poles on 7. Then we can write f(z) = r(z)e ie ^ for some functions r: C — > R and 
9: C — » R, where < r(z) < 00 for all 267, For each z £ C, 9{z) is only defined up to an integer 
multiple of 2ir. However, we can require 6(z) to be continuous as we vary z 6 7, in which case 9{z) 
is determined up to an overall offset by an integer multiple of 2ir. Then for f{z) to be single- valued, 
the function 9{z) must change by some integer multiple of 2tt as z winds around 7. This integer is 
the winding number of f around 7, 



We can compute w^{f) by appealing to the Cauchy argument principle, which says that the 
winding number of / around 7 is determined by the zeros and poles of / inside 7. In particular, 
we have the following: 

Theorem 1 (Cauchy argument principle). Let 7 be a simple, closed, positively-oriented curve in 
C, and let f : C — > C be a meromorphic function with no zeros or poles on 7. Then 




(19) 




(20) 



(21) 



W M) = Z int( 7 )(/) - ^int( 7 )(/) 



(22) 



where Z; nt ( 7 )(/) is the number of zeros of f inside 7, counted with their multiplicity, and -Pmt(7)(/) 
is the number of poles of f inside 7, counted with their order. 
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Levinson's theorem relates the winding of the reflection coefficient around the unit circle T to 
the number of bound states. In particular, we have the following: 

Theorem 2 (Levinson's theorem for graphs). Let R(z) be the reflection coefficient of a graph 
scattering problem with m internal vertices in which there are n c confined bound states, half- 
bound states, and n& evanescent bound states. Furthermore, suppose that either a 7^ or ||6|| 7^ 1 
(or both). Then 

wr(R) = 2(m - n b - n c - \n h ). (23) 

The technical condition on a and b poses no essential difficulty. If a = and ||6|| = 1, then by 
a suitable change of basis, the problem is equivalent to one in which a = and 6 f = (1, 0, . . . , 0). 
In this case, we can effectively view the first vertex of the finite graph as part of the semi-infinite 
path, and we can replace the {m + I) x (m + 1) matrix (4) by the m x m matrix D. Repeating this 
process, we eventually arrive at a graph with a / or ||6|| 7= 1 (or else we arrive at a graph with a 
single vertex, in which case the bound states are easily characterized). 

The proof of Theorem 2 uses the argument principle. Note that since |-R(z)| = 1 for \z\ = 1, 
R{z) has no zeros or poles on the unit circle. Our goal is to count the zeros and poles of R(z) inside 
the unit circle. 

It will be helpful to view R(z) as a quotient as in (12) and consider the numerator and denom- 
inator of this expression separately. Indeed, except at the point z = 0, every zero of R(z) inside T 
corresponds to a zero of Q(z~ 1 ), and every pole of R(z) inside T corresponds to a zero of Q(z). 

Lemma 3. Let < \z\ < 1. Then R(z) = if and only if Q(z~ 1 ) = 0, and R(z)~ 1 = if and only 
ifQ(z) = 0. 

Proof. Suppose that R(z) = 0. We claim that Q(z _1 ) = 0; assume for a contradiction that this is 
not the case. To satisfy R{z) = 0, Q(z) must be unbounded. Since a is a constant and \z\ < 1, 
this can only happen if C(z) is unbounded. However, if C(z) is unbounded then R(z) = —z~ 2 7^ 0. 
This is a contradiction, so we must have Q(z _1 ) = 0. 

Similarly, if i?(z) -1 = and Q(z) 7^ 0, then Q(z~ 1 ) must be unbounded, and since a and z~ l 
are finite, C(z~ 1 ) = C(z) must again be unbounded, showing that R(z)~ l = —z 2 7^ 0. This is a 
contradiction, showing that Q(z) = 0. 

Clearly, if Q{z~ l ) = and Q(z) ± then R(z) = 0, and if Q{z) = and Qiz" 1 ) ^ 
then Riz)" 1 = 0. It remains to consider the possibility that Q(z) = Q(z x ) = 0. In this case, 
z{a + C{z)) = z~ l {a + C(z" 1 )) = z^{a + C{z)). Since a + C(z) + (as a + C(z) = implies 
Q(z) = Q(z^ 1 ) = 1 7^ 0), this shows that z 2 = 1, which violates the condition \z\ < 1. □ 

Thus it suffices for us to count the zeros of Q(z _1 ) and of Q(z) inside T. First we count the 
zeros of Q(z) inside T. We saw in Section 2.2 that every bound state corresponds to a root of Q(z). 
In particular, for any z = ±e~ K € [—1,0) U (0, 1] with Q(z) = 0, we can solve (17) to find a bound 
state. Next we show that zeros of Q(z) on and inside the unit circle can only occur on the real 
axis, which demonstrates that every zero of Q(z) on and inside the unit circle (except at z = 0) 
corresponds to a bound state. 

Lemma 4. Let \z\ < 1 and Q(z) = 0. Then lm(z) = 0. 

Proof. The claim is trivial for z = 0, so we can assume \z\ > 0. Then 

= hn.{z~ x Q{z)) (24) 
= Im(z" 1 ) -Im(C(z)). (25) 
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With z = re %6 we have Im(z 1 ) = —r 1 s'm9. Let D = ^\ \P\ be the spectral decomposition of 
D, where P\ projects onto the eigenspace of D with eigenvalue A. Then we have 

A 



Therefore, 

lm(C(z)) = - V HiWV-Qsiiifr 

1 1 ;J [( r + r -i) cos0 _A]2 + [( r _ r -i) s ine] 2 ' 



Now assume for a contradiction that sin# ^ 0. Then we have 

r -i = (r _ r -i)Y^ \\ p x b \\ 2 (28) 

{ [(r + r- 1 ) cos e- A] 2 + [{r-r~ 1 ) sin 6>] 2 ' 1 j 

A 

For < r < 1, the left hand side is positive but the right hand side is non-positive. This is a 
contradiction, so we must have sin# = 0, i.e., Im(z) = 0. □ 

Furthermore, each bound state contributes a zero of multiplicity one. 

Lemma 5. Suppose that Q(x) = where < \x\ < 1. Then 4~Q(x) ^ 0. 

Proof. By Lemma 4, x G R. Since = 0, a + C(x) = x -1 . Then we have 

^-Q(x) = -(a + C(x))-x^-C(x) (29) 
dx ax 

' S J«SL (30) 
dx ^ x + x" 1 - A 

A 



IPa^H 2 ^- 1 



" x 1 E "( x A +x-i - A) 2 ' ' (31) 

A 

which is strictly negative for < x < 1 and strictly positive for — 1 < x < 0. □ 

So far, we have ignored the point z = 0. Although there are not bound states at z = 0, i?(z) 
can nevertheless have a pole at z = 0. 

Lemma 6. If a ^ 0, £/ien -R(-z) /ias a simple pole at z = 0. If a = and \\b\\ ^ 1, i/ien -R(z) /ias 
no zeros or poles at z = 0. 

Proo/. We have = -Q(z" 1 )/Q(z). Note that C(0) = 0. Since Q(0) = 1, Q(ar) contributes no 
zeros or poles at z = 0. Now Q(z _1 ) = z~ l (z — a — C(z)). If a 7^ 0, then lim z ^oz — a — C(z) = 
—a 7^ 0, so R(z) has a simple pole at z = as claimed. On the other hand, if a = 0, then 
lim 2 _5.o Q(-2 _1 ) = lim z _>o(l — z C(z)) = 1 — \\b\\ 2 . Thus if a = and \\b\\ 7^ 1, has no zeros or 
poles at z = as claimed. □ 

We must also count the zeros of Q(z~ 1 ) inside T. We do this by relating the zeros of Q(z _1 ) 
inside T to the zeros of Q(z) outside Y. Observe that if Q(z^ n 1 ) = where \z- m \ < 1, then Q(z out ) = 
for some z ou t with \z out \ > 1, namely z ou t = z^~. Thus, it will be sufficient to determine the total 
number of zeros of Q{z). 

To count the zeros of Q(z), we use the fact that it is a rational function of z and determine the 
degree of its numerator. First we consider the case where a / 0. 
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Lemma 7. If a ^ 0, then Q{z) is a rational function of z whose numerator has degree 2m — 2n c + l. 
If a = and \\b\\ ^ 1, then Q(z) is a rational function of z whose numerator has degree 2m — 2n c . 

Proof. Using (26), we have 

Q(z) = l-az-z>Y,JP£^- 1 (32) 

A 

(1 - az) I\x(z 2 -Xz + l)-z 2 EaII^II 2 ELyA^ ~ A'z + 1) 

nA(- 2 -A, + i) • (33) 

If A = ±2, then z 2 — Xz + 1 = (z =F l) 2 has two identical roots (z = ±1). For A ^ ±2, we claim 
that the two roots of z 2 — Xz + 1 are distinct. Furthermore, for A ^ A', z 2 — Xz + 1 and z 2 — X'z + 1 
have no common roots. To see this, write z 2 — Xz + 1 = (z — ^+){z — //_) where 



»±:=$±y/{$) 2 -l. (34) 

Observe that if |A| < 2 then fi + is on the half of the unit circle in the complex plane with positive 
imaginary part and /i_ is on the half of the unit circle with negative imaginary part. If A > 2 
then E (1, oo) and fi- € (0,1), and if A < — 2 then fi + € (—1,0) and |U_ € ( — oo, — 1). Thus, 
clearly fj, + ^ fi'_ regardless of whether A = A'. We also claim that ^ fj,'_±_ and /i„ ^ jj!_ when 
A A'. For |A| < 2 this is obvious, simply by comparing real parts. For |A| > 2, observe that /i + 
is a strictly increasing function of A for A > 2 and a strictly decreasing function for A < — 2, and 

is a strictly decreasing function of A for A > 2 and a strictly increasing function for A < — 2; 
therefore, for A ^ A', fj, + ^ fi' + and fj,- ^ fj,'_- 

Now let J(z) denote the numerator of (33), so that 

°< 2 > - rnZ-l+D - (35) 

If a ^ 0, then J{z) is a polynomial with leading term — az 2m+1 , where m is the number of distinct 
eigenvalues of D. On the other hand, if a = and ||6|| ^ 1, then the leading term of J(z) is 

(±-Ex\\Pxb\\ 2 )z 2rh = (l-\\b\\ 2 )z 2rh - 

It remains to count common factors between the numerator and denominator of (33). For an 

eigenvalue A ^ ±2 of D, let fi ^ ±1 be either of the two roots of z 2 — Xz + 1. Then [i is a simple 

root of IIa'^ 2 — A'z + 1), and 

J(/i) = -AWt H (f, 2 - A'/i + 1), (36) 
ava 

so fj, is a root of J(z) if and only if ||-P\&| =0. If A = ±2 is an eigenvalue of D, then fi = ±1 is a 
double root of IIa'O 2 ~ X ' z + If W P M + then ( 36 ) shows that 

\x is not a root of J(z), and 

if ||Pa&|I = then it is easy to see that fi is a root of J(z) of multiplicity at least 2. Thus we see 
that the number of common factors between the numerator and denominator of (33) is twice the 
number of values A such that ||P\6|| = 0. Since the number of confined bound states corresponding 
to A is either the multiplicity of A (if ||Pa&|| = 0) or the multiplicity of A minus 1 (if H-Pa^H 7^ 0), 
the total number of common factors is 2(n c + m — m). Then the degree of the numerator of Q(z) 
is 2m — 2n c + 1 if a ^ and 2m — 2n c if a = and ||6|| ^ 1, as claimed. □ 

Finally, we are ready to prove the main result. 
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Proof of Theorem 2. For any S C C and any meromorphic function /, let Zg(f) denote the number 
of zeros of / that belong to S (counted with their multiplicity) and let Ps(f) denote the number 
of poles of / that belong to S (counted with their order). Let int(r) denote the interior of T, and 
let ext(r) denote the exterior of T (both open sets, not including points on T). Since R(z) is a 
rational function of z, in particular it is a meromorphic function of z; since |-R(z)l = 1 for \z\ = 1, 
it has no zeros or poles on V. Thus, by the argument principle (Theorem 1), 

w T (R) = Z int{r) (R) - P mt[T) (R). (37) 

By Lemma 6, R{z) has a simple pole at z = if a ^ and no zeros or poles at z = otherwise, so 

w r (R) = Z int(r)m (R) - P inm \ { o}(R) ~ S[a + 0] (38) 

where S[a ^ 0] = 1 if a ^ and 5[a ^ 0] = if a = 0. 
By Lemma 3, 

z mt(r)\{o}(R) = ^int(r)\{o}(Q(^ -1 )) (39) 
and P in t(r)\{o}(R) = Z- mt{r )\ {0} (Q(z)). (40) 

In particular, 

w r (R) = ^ nt (r)\{o}(Q(^ 1 )) - Z inm \ {0} (Q(z)) - S[a + 0]. (41) 

By Lemma 4 and Lemma 5, Z int ^\^ y(Q(z)) = Z int ^(Q(z)) = n&. Also, as discussed above, 
we can relate zeros of Q(z~ 1 ) inside F to zeros of Q(z) outside F. In particular, 

^in t( r)\{o}(Q(^" 1 )) = Z cxt(r) (Q(z)) (42) 
= Z c (Q(z)) - Z T (Q(z)) - Z int(r) (Q(z)). (43) 

To compute Zr(Q(z)), recall that zeros of Q(z) at z = ±1 correspond to half-bound states. By 
Lemma 4, these are the only zeros with \z\ = 1, and by Lemma 5, they occur with multiplicity one. 
Thus we have Zr(Q(z)) = n^, the total number of half-bound states. 

Finally, by Lemma 7, Zc(Q(z)) = 2m — 2n c + 5[a ^ 0]. Combining these expressions, we have 

wt(R) = 2m - 2n b - 2n c - n h (44) 
as claimed. □ 



4 Discussion 

We conclude by discussing several possible directions for future work. 

First, it might be interesting to investigate possible applications of bound states to quantum 
computation. Existing quantum algorithms based on graph scattering [2, 3] rely solely on the 
properties of scattering states, but one can imagine that bound states could also encode information 
about the solution of a computational problem. While it might be possible to extract information 
about bound states directly using phase estimation, it could also be useful to infer such information 
indirectly via Levinson's theorem. 

One could also consider generalizing Theorem 2 to the case where n > 1 semi-infinite paths 
are attached to n vertices of an (n + m)-vertex graph. In general, the scattering is described by 
an n x n matrix, the 5-matrix, instead of a single reflection coefficient (see Ref. [2] for details). 
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Similar calculations to those in Section 2.1 can be used to determine the S- matrix in this more 
general case: for a weighted adjacency matrix of the form where A E C nxn , B € C mxn , 

and D £ c mxm ) the S-matrix at momentum fe is S(e lk ), where S(z) = —Q(z~ 1 )/Q(z) with Q(z) = 
I — z(A + B^[(z + z~ 1 )I — D]~ 1 B). Furthermore, similar conditions to those described in Section 2.2 
can be used to characterize the bound states. We expect the number of bound states in such a 
scattering problem to be related to the winding number of det S(z), but we leave the details as a 
topic for future work. 

Finally, the original theorem of Levinson for a potential on a half-line not only describes a 
relationship between the number of bound states and the phase of the reflection coefficient, but also 
sheds light on the extent to which the potential can be inferred from the behavior of the reflection 
coefficient. It might be fruitful to consider such inverse problems in the setting of scattering on 
graphs. The weighted graph corresponding to a given form of the reflection coefficient is certainly 
far from unique, since one may perform arbitrary similarity transformations on the m internal 
vertices without affecting the scattering behavior, but one might attempt to characterize the graph 
corresponding to a given reflection coefficient up to such changes of basis. 
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